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Abstract 

Background: Rice is one of the major crop species in the world helping to sustain approxinnately half of the global 
population's diet especially in Asia. However, due to the innpact of extrenne climate change and global warming, 
rice crop production and yields may be adversely affected resulting in a world food crisis. Researchers have been 
keen to understand the effects of drought, temperature and other environmental stress factors on rice plant growth 
and development. Gene expression microarray technology represents a key strategy for the identification of genes 
and their associated expression patterns in response to stress. Here, we report on the development of the rice 
OneArray® microarray platform which is suitable for two major rice subspecies, yapon/ca and indica. 

Results: The rice OneArray® 60-mer, oligonucleotide microarray consists of a total of 21,179 probes covering 20,806 
genes of joponica and 13,683 genes of indica. Through a validation study, total RNA isolated from rice shoots and 
roots were used for comparison of gene expression profiles via microarray examination. The results were submitted 
to NCBI's Gene Expression Omnibus (GEO). Data can be found under the GEO accession number GSE50844 
(http://www.ncbi. nlm.nih.gov/geo/query/acc.cgi?acc=GSE50844). A list of significantly differentially expressed genes 
was generated; 438 shoot-specific genes were identified among 3,138 up-regulated genes, and 463 root-specific 
genes were found among 3,845 down-regulated genes. GO enrichment analysis demonstrates these results are in 
agreement with the known physiological processes of the different organs/tissues. Furthermore, qRT-PCR validation 
was performed on 66 genes, and found to significantly correlate with the microarray results (R = 0.95, p < 0.001 ^^^). 

Conclusion: The rice OneArray® 22 K microarray, the first rice microarray, covering both japonica and indica 
subspecies was designed and validated in a comprehensive study of gene expression in rice tissues. The rice 
OneArray® microarray platform revealed high specificity and sensitivity. Additional information for the rice 
OneArray® microarray can be found at http://www.phalanx.com.tw/index.php. 
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Background will have to produce 70% more food by 2050 to feed a 

Rice is one of the most important crops in the world - a projected extra 2.3 billion people. As such, rice crop 

staple food supporting more than half of the world s 7 production will play an important role to maintain food 

billion people. By 2050, the global population is antici- security in the coming future. 

pated to expand between 7.5 and 10.5 billion with the In recent years due to abnormal climate changes, sev- 

growth concentrated mainly in rice consuming coun- eral grain countries such as Australia, Brazil, and 

tries. According to a 2009 report by the United Nations Thailand, have suffered frequently from devastating 

Food and Agriculture Organization (FAQ), the world floods and droughts, resulting in global grain crop losses 

and food price inflation. Additionally, extreme weather 
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the Asia Pacific region; for example, Typhoon Morakot 



Hsinchu 30078, Taiwan, Republic of China brought catastrophic damage to Taiwan in 2009 [1]. 
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Furthermore, substantial decreases in rice yields due to 
increased nighttime temperatures associated with global 
warming [2], and increased minimum air temperatures 
during growing seasons have been reported in China and 
Philippines with predictions of this phenomenon continu- 
ing [2,3]. Increasing in atmospheric brown clouds and 
greenhouse gas has been proposed to reduce historical rice 
harvests in India well below expected levels. These studies 
have profound implications for ongoing and future efforts 
for climate and air quality improvements [4,5]. 

Therefore, the development of new rice strains against 
the threat of climate change and water shortages is an 
important issue of food security for the coming future 
[6,7]. Traditional rice breeding techniques required up- 
wards of 10 years to develop new rice strains due to a 
multi-generational process of selecting and preserving 
strain variety in the offspring. The application of mo- 
lecular marker-assisted breeding will facilitate the pyra- 
miding of desirable alleles at multiple loci and shorten 
the time needed for developing new varieties [8]. An in- 
crease in the fundamental knowledge of rice biology (e.g. 
seed formation, disease resistance, growth etc.) is neces- 
sary to provide for effective strategies to improve crop 
yield and production. Nevertheless, the majority of the rice 
consumption is produced by indica subspecies, but the 
greater part of the genomic work is done on japonica. To 
overcome the imbalance in rice related studies, it is im- 
portant to generate a microarray using in both subspecies. 

To aid rice researchers and plant biologists, the rice 
OneArray® gene expression microarray was developed by 
Phalanx Biotech Group Inc. Based on the rice genome 
sequences from the Rice Genome Annotation Project 
(version 6.1) and the Beijing Genome Institute 2008 
database, probes were designed and selected to cover 
90% of the well-annotated Gene found on both japonica 
and indica subspecies. Researchers will find the rice 
OneArray® microarray suitable for large-scale basic stud- 
ies, stress physiology research, and biomarker discovery. 

In this paper, we report the development and validation 
of the 22 k rice OneArray® oligo- microarray platform 
manufactured by Phalanx Biotech Group, Inc. The high 
data reproducibility on technical replicates and the plat- 
forms capacity to identify differential gene expression in 
different tissues and rice subspecies was demonstrated. 

Results and discussion 

Rice OneArray® gene selection and probe design 

The primary goal of this project was to develop a rice 
microarray platform to study gene expression patterns 
relevant to important biological and physical controls 
across japonica and indica subspecies. The rice OneArray® 
microarray was specifically designed to cover important 
regulatory pathways and to include genes involved in the 
biological function of chloroplasts, oxidative stress, grain 



quality, nitrogen phosphate, sugar synthesis, photosyn- 
thesis, plant hormone, anther development, and transcrip- 
tion factors. 9719 japonica target genes were curated and 
selected from a compilation of databases including 
Michigan State University database (ftp://ftp.plantbiology. 
msu.edu/pub/data/Eukaryotic_Projects/o_sativa/annotation_ 
dbs/pseudomolecules/chloroplast.dir/chrC.cDNA), GOSlim 
(http://rice.plantbiology.msu.edu/annotation_pseudo_ 
goslim.shtml), K E G G (http://www.kegg.jp/kegg-bin/ 
show_organism?menu_type=pathway_maps&org=osa), 
Gramene (ftp://ftp.gramene.org/pub/gramene/pathways/ 
ricecyc/), and PlnTFDB (http://plntfdb.bio.uni-potsdam. 
de/v3.0/). mRNA transcript sequences of each japonica 
target gene were first subjected for microarray probe de- 
sign by using IMPORT software (Industrial Technology 
Research Institute of Taiwan, R.O.C). Probes were de- 
signed according to the following criteria; 60 nucleotides 
in length, GC% between 40-60%, fewer than 6 simple nu- 
cleotide repeats, and probe location within 1200 bp from 
3' terminus. To remove the probes with non-specific bind- 
ing and strong secondary-structure, probe sequences were 
run through Blast analyses against japonica (from Rice 
Genome Annotation Project (version 6.1)) and indica 
whole genome sequences (from BGI; Beijing Genome In- 
stitute). Based on the above probe design procedure, 92% 
of 9719 japonica target genes and 86% of 8995 indica 
target genes were selected for inclusion on the rice 
OneArray® microarray. The probe set was further supple- 
mented with probes against well- annotated genes as de- 
fined by Gene Ontology. In total, 21,179 probes (Figure 1) 




Figure 1 Number of probes In the microarray represented in 

japonica and indica subspecies. A total of 21,179 probes were 

classified into seven groups (A, B, C E, F and G) based on Blast 
analysis against corresponding genomes of both subspecies. For 
japonica, 20,806 probes were represented among A, B, D, E, F, and G 
groups, with single hit probes in group A, B, and D, and multiple-hits 
probes in group E, F, and G. For indica, 13,683 probes were incorporated 
in A, C, E, D, and G groups, with single hit probes in group A, C, and E 
and multiple hits probes in group D and G group. 

V ) 
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were selected and designed based on the rice whole gen- 
ome sequences from japonica and indica subspecies, plus 
824 control probes including IHC (Intrinsic Hybridization 
Control), IHL (Intrinsic Hybridization Ladder), ITQC (In- 
trinsic Target Quality control), and negative controls. In 
summary, the rice OneArray® microarray is suitable for the 
detection of 20,806 genes of japonica and 13,683 genes 
of indica. 

Array quality 

Each microarray undergoes a spot QC process to evaluate 
probe deposition and immobilization efficiency. In brief, 
microarrays were incubated with random, 10-mer oligo 
probes labeled with Cy3 using the standard hybridization 
protocol (Figure 2A). The sensitivity and dynamic detec- 
tion range of the rice OneArray miroarray were tested 
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Figure 2 Validation of Rice OneArray ® microarray quality. A. 

Examination of spot quality control. B. Sensitivity testing by using 
five external spike-in probes. 



using commercially available external control probe sets 
(external spike-in system) supplemented with 10 (ig Cy5- 
labeled aRNA of rice shoot. The results showed that the 
minimal detectable concentration of probes was approxi- 
mately 0.05 pM (Figure 2B). In summary, the rice micro- 
array demonstrated high sensitivity and dynamic range in 
this gene expression profiling study. 

Rice OneArray® microarray technical performance 

The Phalanx microarray platform is based upon the 
hybridization of a single labeled sample (derived from 
RNA), followed by one-channel detection. The intensity 
of the hybridization signal is used to determine target 
concentration. In order to validate the technical quality 
of each probe in our arrays, we carried out 10 independ- 
ent hybridizations on samples representing two different 
rice tissues - root and shoot. 

To examine the gene expression profiles between rice 
root and shoot development, total RNA extracted 
from rice root and shoot were processed on the rice 
OneArray® microarray following the standard protocol 
using five arrays for each tissue type. Raw expression 
data from 10 microarrays (e.g. 5 arrays x 2 tissues) were 
normalized and Pearsons correlation coefficients were 
calculated for the data sets of hybridization signal inten- 
sities. All normalized and raw data were submitted to 
NCBIs Gene Expression Omnibus (GEO) for others to 
examine. The data are accessible via GEO Series acces- 
sion number GSE50844 http://www.ncbi.nlm.nih.gov/ 
geo/query/acc.cgi?acc=GSE50844. It was demonstrated 
that the average spot number of each tissue was approxi- 
mately 18,000 spots with an average signal intensity of 
each spot of 3,000. High correlation coefficients were ob- 
tained in all cases and the results obtained for the rice root 
and shoot sample are shown in Figure 3A (r = 0.996; 
p-value <0.001***) and 3B (r = 0.998; p-value <0.001***) re- 
spectively. Furthermore, the significant correlation was ob- 
served between technical repeats (Table 1, R> 0.983) in 
each of the 5 arrays. In summary, 90% of probes can be 
detected in rice root and shoot tissue, the detectable spot 
percentage is higher than other species array data (60- 
70%). The higher coverage rate of the rice OneArray® 
microarray may be attributed to the comprehensive target 
gene selection of genes relevant to rice development. 
These results also demonstrate high correlation between 
different technical experiments underlining the high preci- 
sion manufacturing of the rice OneArray® microarray 
platform. 

Comparison of gene expression in rice root and shoot 
using gene ontology analysis 

To elucidate the genes regulating rice tissue develop- 
ment, comparisons between shoot and root gene expres- 
sion profiles from the rice OneArray® microarray were 
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Figure 3 Rice OneArray® microarray reproducibility. Scatter plot 
indicating correlation between microarray technical replicates. High 
correlations were found between all technical replicates. A. rice 
shoot (r = 0.996, p < 0.00 r**). B. rice root (r = 0.998, p < 0.001 



normalized and analyzed using Rosetta software. A list 
of differentially expressed genes was generated; among 
3,138 up-regulated genes (Additional file 1), 438 were 
shoot-specific genes (Additional file 2), and among 3,845 
down-regulated genes (Additional file 3), 463 were root- 
specific genes (Additional file 4). Gene set analysis was 
performed using Gene Ontology terms with functional 
annotation, as described in DAVID Bioinformatics Re- 
sources 6.7 (http://david.abcc.ncifcrf.gov/) [9,10]. First, 
this method detects significantly up- or down- regulated 
clusters of functionally related genes in lists ordered by 
differential expression. Annotated genes in these differ- 
ent groups were then classified into different GO bio- 
logical processes and the percentages of differential gene 
expressions were calculated for each process. Among 
up-regulated genes, GO biological processes included 
oxidation reduction (4.33%), photosynthesis (1.15%), 
pigment metabolic (0.41%) and biosynthetic process 
(0.38%), and fatty acid metabolic (0.48%) and biosyn- 
thetic processes (0.45%). A similar set of GO processes 
were observed in the shoot-specific cluster (Table 2). 
Among down-regulated genes, GO biological processes 
included regulation of transcription (2.86%), response to 
oxidative stress (0.7%), and cell wall polysaccharide 
metabolic processes (0.13%), of which the majority was 
observed in root-specific cluster (Table 3). Overall, these 
results are in general agreement with the known physio- 
logical processes of the different organs/tissues suggest- 
ing the rice OneArray® platform is capable of providing 
reliable gene-expression data. 

qRT-PCR validation 

To further validate the microarray results, quantitative real 
time PGR (qRT-PCR) assays were performed on the same 
RNA samples used for microarray analysis. A total of 66 
genes at varying expression levels including up-regulated, 



Table 1 The correlation between technique repeats of rice shoot and root 



Shoot Root 







HOI H02 


H03 


H04 


H05 


H06 


H07 


H08 


H09 


H10 




HOI 


1.000 0.989 


0.987 


0.986 


0.985 


0.670 


0.673 


0.674 


0.675 


0.672 




H02 


1.000 


0.991 


0.985 


0.988 


0.678 


0.678 


0.681 


0.682 


0.678 


Shoot 


H03 




1.000 


0.987 


0.988 


0.681 


0.679 


0.683 


0.681 


0.678 




H04 






1.000 


0.983 


0.667 


0.673 


0.671 


0.675 


0.671 




H05 








1.000 


0.674 


0.679 


0.677 


0.681 


0.677 




H06 










1.000 


0.984 


0.994 


0.987 


0.983 




H07 












1.000 


0.987 


0.993 


0.985 


Root 


H08 
H09 
H10 














1.000 


0.990 
1.000 


0.985 
0.986 
1.000 
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Table 2 Go term significantly represented in up-regulated and shoot-specific gene clusters 





Term (biological process) 


Count 


Percentage 


P-value 


Up-regulated genes: 3138 (shoot vs. root) 


Oxidation reduction 


136 


4.33% 


l.lOE-17 




Photosynthesis 


36 


1.15% 


1.56E-15 




Photosynthesis, light harvesting 


10 


0.32% 


2.70E-07 




Pigment metabolic process 


13 


0.41% 


6.27E-06 




Pigment biosynthetic process 


12 


0.38% 


9.67E-06 




Photosynthesis, light reaction 


12 


0.38% 


1 .30E-05 




Generation of precursor metabolites and energy 


36 


1.15% 


5.09E-05 




Fatty acid biosynthetic process 


14 


0.45% 


0.00286 




Fatty acid metabolic process 


15 


0.48% 


0.004194 




Electron transport chain 


14 


0.45% 


0.006286 




Aminoglycan metabolic process 


7 


0.22% 


0.006382 




Lipid biosynthetic process 


21 


0.67% 


0.007102 




Cellular aldehyde metabolic process 


4 


0.13% 


0.007622 




Lipid transport 


12 


0.38% 


0.009804 




Lipid localization 


12 


0.38% 


0.009804 


Shoot-specific genes: 438 


Oxidation reduction 


40 


9.13% 


3.40E-07 




Fatty acid metabolic process 


8 


1 .83% 


1 .30E-03 




Recognition of pollen 


6 


1 .37% 


2.70E-03 




Pollination 


6 


1 .37% 


2.70E-03 




Cell recognition 


6 


1 .37% 


2.70E-03 




Pollen-pistil interaction 


6 


1 .37% 


2.70E-03 




Fatty acid biosynthetic process 


7 


1 .60% 


2.80E-03 




Reproductive cellular process 


6 


1 .37% 


2.90E-03 




Regulation of transcription 


30 


6.85% 


l.lOE-02 



down-regulated, and not differentially-expressed were se- 
lected for validation by comparison of root and shoot ex- 
pression profiles. Comparisons between microarray and 
qRT-PCR data are shown (Table 4), and microarray results 
correlated well with qRT-PCR validation (Figure 4, R = 
0.95, p < 0.001**^^). 

Conclusions 

A newly-designed rice microarray, the rice OneArray® 
22 K microarray, was provided for examining both 
japonica and indica subspecies. It was demonstrated this 
platform displayed high specificity and sensitivity follow- 
ing a comprehensive validation. Based on the unique de- 
sign, we believe this microarray will be of interest to many 
researchers in rice studies, especially in important bio- 
logical and physical controls, and it can be used to facili- 
tate the functional studies toward a hybrid subspecies. 

Methods 

Tissue preparation and total RNA extraction 

A three-leaf-stage japonica subspecies (Tainung 67, 
TNG 67) was selected and subjected for total RNA 



extraction. In general, 100 mg of rice tissue was cut into 
5 cm lengths and stored immediately in RNAlater (Invi- 
trogen, Carlsbad, CA, USA) at 4°C until RNA isolation. 
Rice tissues were homogenized using a RNase-free mor- 
tar before performing RNA extraction, and total RNA 
was isolated from rice roots and shoots using the Qiagen 
RNeasy Mini kit (Qiagen, Chatsworth, CA, USA) ac- 
cording to manufacture s protocols. 

cRNA amplification 

1 (ig of total RNA was converted to double stranded 
cDNA using reverse transcriptase, and amplified by 
in vitro transcription using MessageAmpII aRNA Amplifi- 
cation kit (Ambion Inc., Austin, Texas, USA). The synthe- 
sized cRNA was subsequently conjugated with Cyanine 5 
NHS ester dye (GE Healthcare, Milwaukee, WI, USA). 
cRNA yield and labeling efficiency was calculated based on 
ND-1000 spectrophotometer measurements (NanoDrop 
Technologies, Wilmington, DE, USA). Incorporation rates 
of 20-60 dye molecules per 1,000 bases (20-33 bases/dye 
molecule) yielded the most usable data. 



Chen et al. BMC Research Notes 2014, 7:15 
http://www.bionnedcentral.conn/1 756-0500/7/1 5 



Page 6 of 9 



Table 3 Go term significantly represented in down-regulated and root-specific gene clusters 





Term (biological process) 


Count 


Percentage 


P-value 


Down-regulated genes: 3845 (shoot vs. root) 


Regulation of transcription 


110 


2.86% 


6.28E-13 




Regulation of transcription, NA-dependent 


71 


1 .85% 


1.20E-10 




Regulation of RNA metabolic process 


71 


1 .85% 


1.48E-10 




Transcription 


65 


1 .69% 


8.87E-10 




Response to oxidative stress 


27 


0.70% 


1.93E-08 




Oxidation reduction 


81 


2.11% 


1.21 E-05 




Cell wall polysaccharide metabolic process 


5 


0.13% 


0.008210995 




Oxylipin metabolic process 


5 


0.13% 


0.008210995 




Xylan metabolic process 


5 


0.13% 


0.008210995 




Xylan catabolic process 


5 


0.13% 


0.008210995 




Oxylipin biosynthetic process 


5 


0.13% 


0.008210995 


Shoot-specific genes: 438 


Response to oxidative stress 


11 


2.51% 


4.30E-05 




Oxidation reduction 


27 


6.16% 


1 .30E-03 




Cell wall polysaccharide metabolic process 


4 


0.91% 


1 .50E-03 




Xylan metabolic process 


4 


0.91% 


1 .50E-03 




Hemicellulose metabolic process 


4 


0.91% 


1 .50E-03 




Xylan catabolic process 


4 


0.91% 


1 .50E-03 




Regulation of transcription 


27 


6.16% 


5.30E-03 




Polysaccharide metabolic process 


8 


1 .83% 


6.90E-03 



Microarray pre-hybridization 

Rice OneArray® microarrays were pre-heated at 60°C for 
10 min in hybridization oven. Microarray slides were 
placed inside a falcon tube containing 100% ethanol, in- 
cubated for approximately 15 sec, shaken for 20 sec, and 
thoroughly rinsed with deionized water to remove any 
residual ethanol Next, the microarray slides were fully 
submerged in an abundant amount of pre-hybridization 
solution (5X SSPE, 0.1% SDS, and 1% BSA) for 1 hr at 42°C. 
After 1 hr, slides were transferred to room-temperature dis- 
tilled water and washed gently for 2 min. Slides were spun 
dry for 2 min and stored in a dry and dark place until 
hybridization. 

Microarray hybridization 

10 (ig of cRNA was fragmented by using RNA Fragmenta- 
tion Reagent kit (AM#8740, Ambion Inc., Austin, Texas, 
USA), and then denatured in a PGR machine at 95°C for 5 
minutes and held at 60°C. Fragmented cRNA was hybrid- 
ized on the rice One Array® (Phalanx Biotech Group, 
Taiwan) at 50°C for 14-16 hrs. After hybridization, the mi- 
croarrays were washed sequentially in 2X SSC containing 
0.2% SDS solution for 5 min at 42°C, 2X SSC for 5 min at 
42°C, and 2X SSC for 5 min at room temperature. Finally, 
the microarrays were spun dry with a centrifuge for at 
least one minute and stored dry in the dark until ready for 
scanning. 



Image scanning 

Raw intensity signals for each scanned microarray were 
captured at 10- (im resolution using GenePix Personal 
4000B (Molecular Devices Corporation, Sunnyvale, CA, 
USA), quantified by GenePix™ Pro 4.0 software (Molecular 
Devices Corporation, Downingtown, PA, USA), and stored 
in GPR format. Microarray images were saved as TIFF 
files. Auto Photomultiplier tube (PMT) settings were se- 
lected and adjusted to include the overall feature inten- 
sities of Cy5 channel. 

Data processing and statistical analysis 

The data from all microarrays was processed using pro- 
prietary modeling techniques developed on the Rosetta 
Resolver® System (Rosetta Biosoftware, Seattle, WA, 
USA). Raw data is comprised of probe intensities, back- 
ground values, detected signals, signal-to-noise ratio 
data, probe identification and gene annotations. After 
probe filtering based on flag note criteria, normalization 
of raw intensity was achieved by median scaling and the 
mean of the technical repeats. The log2 (Ratio) were cal- 
culated by pair-wise combination and error weighted 
average. Significant differentially expressed genes (DE 
genes) were selected according to its log2 (Ratio) and 
P-value based on the following criteria; log2 (Ratio) > = 1 
and P-value (differentially expressed) <0.05. 



Chen et at. BMC Research Notes 2014, 7:15 
httpy/www.biomedcentral.com/l 756-0500/7/1 5 



Page 7 of 9 



Table 4 Comparison of expression levels (log2 ration) from qRT-PCR and microarray for selected target genes 



Target gene 



MSU Gene ID 



Microarray value 
(log2ratio) 



qRT-PCR value 
(log2ratio) 



Zinc finger, C3HC4 type, domain containing protein, expressed Os07g29600 0.73 

bZIP transcription factor domain containing protein, expressed Os07g48180 -0.04 

Anl<yrin repeat-ricli protein, putative, expressed Os08g 15840 -0.99 

bZIP transcription factor domain containing protein, expressed Os08g26880 1.42 

Zinc finger family protein, putative, expressed Os08g03310 4.58 

MYB family transcription factor, putative, expressed Os08g33660 3.93 

B-box zinc finger family protein, putative, expressed Os09g35880 0.93 

Zinc finger, C3HC4 type domain containing protein, expressed Os09g33670 -0.13 

MYB family transcription factor, putative, expressed Os09g36250 0.75 

OsWRKY80 - Superfamily of TPs having WRKY and zinc finger Os09g30400 -0.17 
domains, expressed 

BHLH transcription factor, putative, expressed Os09g32510 2.66 

Zinc finger, C3HC4 type domain containing protein, expressed Os09g26400 -0.17 

Dof zinc finger domain containing protein, putative, expressed Os09g29960 2.57 

Auxin response factor 18, putative, expressed Osl0g33940 -1.1 1 

Transcription factor, putative, expressed Osl2gl3170 -0.08 

ZOS12-05 - C2H2 zinc finger protein, expressed Osl2g31840 -0.45 

OsRR4 type-A response regulator, expressed Os01g72330 2.53 

Inducer of CBP expression 2, putative, expressed Os01g70310 4.86 

PHD finger protein, putative, expressed Os01g66420 0.02 

SNF2 family N-terminal domain containing protein, expressed Os02g02290 0.14 

MYB family transcription factor, putative, expressed Os01g59660 -1.41 

No apical meristem protein, putative, expressed Os01g48130 0.31 

GATA transcription factor 25, putative, expressed Os02g05510 -0.53 

PHD finger protein, putative, expressed Os02g35600 2.08 

AP2 domain containing protein, expressed Os02g29550 0.42 

Two-component response regulator, putative, expressed Os02g55320 0.84 

OsMADS22 - MADS-box family gene with MIKCc type-box, expressed Os02g52340 2.78 

Histidine kinase, putative, expressed Os02g50480 0.12 

DIRP family protein, putative, expressed Os03g43800 -0.18 

OsWRKY55 - Superfamily of TPs having WRKY and zinc finger Os03g20550 -2.1 1 
domains, expressed 

ZOS3-23 - C2H2 zinc finger protein, expressed Os03g61 640 -0.64 

NAC domain-containing protein 67, putative, expressed Os03g60080 -4.04 

MYB family transcription factor, putative, expressed Os04g49450 -2.27 

Auxin response factor, putative, expressed Os04g57610 1.17 

MYB_AI protein, putative, expressed Os04g58020 0.1 

Zinc RING finger protein, putative, expressed Os06g03580 -2.82 

bHelix-loop-helix transcription factor, putative, expressed Os05g46370 -0.55 

Transcription factor, putative, expressed Os05g37170 -2.58 

myb-like DNA-binding domain containing protein, expressed Os06g24070 4.09 

Zinc finger protein, putative, expressed Os06g33810 0.27 

N-terminal asparagine amidohydrolase, putative, expressed Os06g41390 -0.6 

TIP41, putative, expressed Os03g55270 -0.73 



0.14 

-0.7 

-0.66 

1.19 

4.99 

2.8 

0.97 

-0.12 

0.2 

-0.48 

1.73 

-0.53 

2.71 

-0.32 

0.66 

-0.38 

2.42 

4.05 

-0.28 

0.11 

-0.76 

-0.02 

-0.09 

0.74 

1.23 

0.46 

2.13 

0.49 

0.21 

-2.45 

0.3 

-4.17 

-2.01 

1.13 

-0.12 

-2.41 

-0.81 

-1.44 

4.63 

-0.25 

-0.02 

-0.5 



Lrien er c//. d/wl nesecircn /vores zu i4, /.id 
httpy/www.biomedcentral.com/l 756-0500/7/1 5 






rage o or y 


Table 4 Comparison of expression levels (log2 ration) from qRT-PCR and microarray for selected target genes 

(Continued) 


B3 DNA binding domain containing protein, expressed 


Os01g04800 


-3.02 


-2.99 


ZOSl-03 - C2H2 zinc finger protein, expressed 


Os01g04120 


-3.41 


-4.35 


TCP family transcription factor, putative, expressed 


Os01gll550 


2.09 


3 


RNA-binding zinc finger protein, putative, expressed 


Os07g48410 


0.68 


0.85 


r^RAQ famik/ trancrripitipin factor Hpimain pontaininpi Pirotoin ovpiroccoH 
vjnAO Idlllliy LI ell loLI ipLIUI 1 IdLLUI UUIIIdlll LUI ILdll III ly piULclll, cApicoocU 


wou/ yoouou 


1 77 


1 44 


Rpcnnnc;p rpniil^tnr rprpi\/pr Hnrn;^in rnnt;^ininn nrntpin PYnrp<^<^pH 


CiqCiJ n4-Q4-f)f) 


-1.54 


-0.69 


AiiYin rpt^nnnQP fpirtnr ni it3ti\/p PYnrpc:cpH 

AAUAIII lt;o|JUIIot; IdLLUI, jJULdLIVc, CAjJI CooCU 


wouoy^u^uu 


—0 82 


0 55 


nui I IcUUUA dooULIdLcU IcULIIIt; Zipptrl, pULdLIVc, t:A|JI coocU 


n<;nQn71 1 RD 


1 06 


0 23 


INU dpiLdI lllcllbLclll piULclll, pULdLIVc, cApicbbcU 


wb 1 uy^z 1 D\j 


— D 
U.Uj 


— D 1 Q 


IVIl D idlllliy LIdl IbLI ipLIUI 1 IdLLUI, pULdLIVc, cApicbbcU 


wbu 1 y / 1 u 


—4 45 




UZ_lr LI dl loLI ipLIUI 1 IdLLUI UUIIIdlll LUI 1 Ldll III ly piULclll, cApi co jfcrU 


n<;ni nssi sn 


0 1 8 


0 1 5 


nc;\A/RKY71 - "^iinprfpimiiv nf TF<; h?i\/inn WRKY ?inrl yinr finnpr 
Wjvvr\i\i/i j)U|-Jciidiiiiiyui iioiidviiiy vvr\i\i diiuz.iiiLiiiiyci 

domains, expressed 


wouzy uoHHU 


—5 14 


—2 75 


Etiiylene-responsive transcription factor, putative, expressed 


Os02g43790 


-4.02 


-1.47 


TCP family transcription factor, putative, expressed 


Os02g42380 


4.84 


4.24 


RING-H2 finger protein, putative, expressed 


Os02g54830 


0.35 


-0.04 


Dof zinc finger domain containing protein, putative, expressed 


Os02g47810 


2.99 


1.97 


Ethylene-responsive transcription factor, putative, expressed 


Os03g09170 


-5.3 


-3.76 


Homeobox associated leucine zipper, putative, expressed 


Os03g08960 


3.17 


2.13 


TCP-domain protein, putative, expressed 


Os02g51280 


0.05 


0.4 


Response regulator receiver domain containing protein, expressed 


Os03g 17570 


-1.86 


-1.42 


ZOS3-21 - C2H2 zinc finger protein, expressed 


Os03g60560 


-5.18 


-6.27 


MYB family transcription factor, putative, expressed 


Os05g04820 


0.68 


0.7 


Actin, putative, expressed 


Os05g36290 


-0.97 


0 


TCP-domain protein, putative, expressed 


Os06g 12230 


0.89 


0.29 



Gene expression was calculated as log2ratio of shoot vs. root samples. Gene ID was obtained from Michigan State University database. 



Quantitative real-time RT-PCR 

Total RNA was isolated according to the methodology pre- 
viously described. Primer designs are listed in Additional 
file 5. 2 ug of total RNA was used to synthesize first-strand 
cDNA using random hexamers and cDNA reverse tran- 
scription kit (Applied Biosystems, Foster City, CA, USA) in 
a reaction volume of 20 [iL, The PGR reactions were per- 
formed in a 20 [iL volume containing Ix Fast SYBR Green 
Master Mix, 20 ng cDNA, 10 nM of forward primers, 10 
nM of reverse primers using an ABI Prism 7900 HT se- 
quence detection system (Applied Biosystems, Foster Gity, 
GA, USA). The amplification conditions were as follows; 
95°G for 20 sec, followed by 40 cycles of 95°G for 5 sec and 
60°G for 30 sec, and final soak at 4°G. Actin RNA served as 
the internal control. The levels of 66 mRNA expressions in 
each of the rice tissue samples were measured by using the 
^-AACt j-Qg^j^Q^ measurements were performed in 
triplicate and the experiments were repeated at least 
twice. 




qRT-PCR data (log2 ratio) 

Figure 4 Correlation between qRT-PCR and Rice OneArray® 
microarray results. Statistically significant correlation (r = 0.947, 
p < 0.001***) was obtained for all 66 tested genes. 
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Additional files 



Additional file 1: Analyses of up-regulated gene expression be- 
tween rice shoot and root tissues. Excel file containing the raw data of 
up-regulated genes represented for shoot vs. root. 

Additional file 2: Shoot-specific expression genes. Excel file 
containing the raw data of shoot-specific expression genes represented 
for comparison between shoot and root. 

Additional file 3: Analyses of down-regulated gene expression be- 
tween rice shoot and root tissues. Excel file containing the raw data of 
down-regulated genes represented for shoot vs. root. 

Additional file 4: Root-specific expression genes. Excel file containing 
the raw data of root-specific expression genes represented for compari- 
son between shoot and root. 

Additional file 5: Primer design list for qRT-PCR validation. Word 
table containing the selected gene description and primer sequences 
indicated the primers used for qRT-PCR examination. 



9. Da Wei Huang BTS, Lempicki RA: Systematic and integrative analysis of 
large gene lists using DAVID bioinformatics resources, Nat Protoc 2008, 
4(1):44-57. 

10. Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward 
the comprehensive functional analysis of large gene lists. Nucleic Acids 
Res 2009, 37(1):1-13. 



doi:1 0.11 86/1 756-0500-7-1 5 

Cite this article as: Chen et al.: Development of a microarray for two 
rice subspecies: characterization and validation of gene expression in 
rice tissues. BMC Research Notes 2014 7:15. 



Competing interests 

The authors declare that they have no competing interests. 



Author's contributions 

SL drafted the manuscript. JSC, SCL, CYC, ^H, PHP, and LC collected the 
genes information and designed those probes. JSC and PHP tested the 
microarray quality and performed microarray experiment for rice tissues. SCL 
and LC performed bioinformatic analysis of microarray QC testing, the 
statistical analysis of expression data and GO functional analysis. JSC and SCL 
designed and coordinated the study. All authors read and approved the final 
manuscript. 



Acknowledgements 

We would like to thank Professor Men-Chi Chang (Department of Agronomy, 
National Taiwan University), Maurice S. B. Ku (Institute of Agricultural Biotech- 
nology, College of Agriculture, National Chiayi University), and Yong-Pei Wu 
(Chiayi Agricultural Experiment Station, Taiwan Agricultural Research Institute) 
for providing rice tissues and advising on the microarray design. This work 
was partially supporting by a grant from the Council of Agriculture, 
Taipei, Taiwan. 



Received: 21 March 2013 Accepted: 3 January 2014 
Published: 8 January 2014 



References 



Butler CD: Food security in the asia-pacific: climate change, phosphorus, 
ozone and other environmental challenges. Asia Pac J Clin Nutr 2009, 
18(4):590-597. 

Peng S, Huang J, Sheehy JE, Laza RC, Visperas RM, Zhong X, Centeno GS, 
Khush GS, Cassman KG: Rice yields decline with higher night temperature 
from global warming. Proc Natl Acad Sci U 5 A 2004, 101(27):9971-9975. 
Tao F, Yokozawa M, Xu Y, Hayashi Y, Zhang Z: Climate changes and trends 
in phenology and yields of field crops in China, 1981-2000. Agr Forest 
Meteorol 2006, 138(l):82-92. 

Auffhammer M, Ramanathan V, Vincent JR: Integrated model shows that 
atmospheric brown clouds and greenhouse gases have reduced rice 
harvests in India. Proc Natl Acad Sci USA 2006, 1 03(52):1 9668-1 9672. 
Cramer W: Air pollution and climate change both reduce Indian rice 
harvests. Proc Natl Acad Sci USA 2006, 1 03(52):1 9609-1 961 0. 
Luo L: Breeding for water-saving and drought-resistance rice (WDR) in 
China. J Exp Bot 201 0, 61 (1 3):3509-35 1 7. 

Redman RS, Kim YO, Woodward CJ, Greer C, Espino L, Doty SL, Rodriguez 
RJ: Increased fitness of rice plants to abiotic stress via habitat adapted 
symbiosis: a strategy for mitigating impacts of climate change. PLoS One 
201 l,6(7):el 4823. 

Collard BCY Vera Cruz CM, McNally KL, Virk PS, Mackill DJ: Rice molecular 
breeding laboratories in the genomics era: current status and future 
considerations. Int J Plant Genonn 2008. doi:l 0.1 155/2008/524847. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



CJ 



BioMed Central 



